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We consider a two-dimensional weakly interacting ultracold Bose gas whose constituents are two- 
level atoms. We study the effects of a synthetic density-dependent gauge held that arises from 
laser-matter coupling in the adiabatic limit with a laser conhguration such that the single-particle 
zeroth-order vector potential corresponds to a constant synthetic magnetic held. We hnd a new 
exotic type of current non-linearity in the Gross-Pitaevskii equation which affects the dynamics of 
the order parameter of the condensate. We investigate the rotational properties of this system in 
the Thomas-Fermi limit, focusing in particular on the physical conditions that make the existence of 
a quantized vortex in the system energetically favourable with respect to the non-rotating solution. 
We point out that two different physical interpretations can be given to this new non linearity: hrstly 
it can be seen as a local modihcation of the mean held coupling constant, whose value depends on 
the angular momentum of the condensate. Secondly, it can be interpreted as a density modulated 
angular velocity given to the cloud. Looking at the problem from both of these viewpoints, we show 
that the effect of the new nonlinearity is to induce a rotation to the condensate, where the transition 
from non-rotating to rotating states depends on the density of the cloud. 

PACS numbers: 


I. INTRODUCTION 

Since its theoretical prediction and first experimental 
realization in 1995 [H, [Ij, Bose-Einstein condensation in 
atomic gases has drawn the attention of the scientific 
community. One of the most intriguing phenomena that 
have been studied since the realization of these macro¬ 
scopic quantum fluids, and in which the signature of 
quantum physics is evident, is represented by the quan¬ 
tized circulation of the velocity field [^-Q. The feature 
that differentiates most a quantum fluid from a classi¬ 
cal one, is in the way it responds to rotation. In con¬ 
trast to a classical fluid, which rotates as a rigid body, 
a Bose-Einstein condensate rotates by forming quantized 
vortices. 

The physics of rotating superfluids and their vortex 
structure is well studied JM, ranging from early studies 
of superfluid Helium Isl r9|| to vortices in strongly inter¬ 
acting quantum gases [I^. Moreover, strongly correlated 
phases in the limit of high vortex numbers, which can be 
regarded as bosonic quantum Hall states, can also ap¬ 
pear in rotating condensates EH- Rotating a condensate 
does however have its limitations. It is only equivalent 
to a uniform magnetic field, and it is not always an easy 
task to control and rotate a small number of particles 
which is needed if one wants to obtain a large flux per 
particle and study quantum Hall physics. For sufficiently 
rapid rotation, or strength of the synthetic magnetic field, 
one hnds that the ground state of the Bose-Einstein con¬ 
densate will correspond to a vortex state. The exper¬ 
imental progress, however, in realising rapidly rotating 
superfluids with ultracold atoms came by taking a differ¬ 
ent route. It was the emulation of orbital magnetism for 
neutral atoms taking advantage of laser-matter coupling 
that made it possible to study magnetic effects without 
physically rotating the gas [1^ [1^ . This has opened up 


a plethora of studies, both in terms of experiments and 
theorw which cover phenomena such as quantised vor¬ 
tices Q, spin-orbit coupling [l^ - [T7j| and various concepts 
from mathematical physics and high energy physics [18l - 
[ 2 ^ . In the standard scenario with synthetic gauge poten¬ 
tials, the orbital magnetism stems from the interaction 
between light and matter. This means that the emerg¬ 
ing gauge potentials are static, and are only given by the 
parameters of the incident lasers. In other words, there 
is no back-action between the gauge held and the mat¬ 
ter held. A non-linear gauge potential on the other 
hand, would change this fact, and result in some intrigu¬ 
ing dynamics [^ . 

In this paper we will address one basic question: how 
are the ground state properties and vortex solutions af¬ 
fected by the presence of a nonlinear gauge potential 
which is proportional to the density of the condensate? 
We will consider an atomic Bose-Einstein condensate 
subject to a laser held conhguration which gives rise to a 
synthetic vector potential, whose zeroth-order contribu¬ 
tion with respect to the ratio between the atom-atom 
and atom-held interaction, emulates a constant mag¬ 
netic held. However, in a mean held picture, detun¬ 
ings from the atomic resonance can be induced by colli¬ 
sions, which results in an unconventional current nonlin¬ 
earity and corresponds to the semiclassical, 

or mean-held limit of an interacting gauge theory [ 2 ^. 

In two dimensions, we show that the current non¬ 
linearity gives rise to a non-trivial non-linear term, pro¬ 
portional to the angular momentum of the condensate. 
It is different from the standard nonlinearity appearing 
in the Gross-Pitaevskii (GP) equation, which affects the 
rotational properties of the condensate in an indirect way 
by modifying the structure of a vortex state, in particu¬ 
lar the size of the vortex core which is of the order of the 
healing length of the condensate, and the overall width 
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of the cloud. We will show that the new nonlinearity 
mentioned above has direct effects on the vorticity of the 
system, because of its explicit dependence on the angu¬ 
lar momentum. We address in particular the problem of 
the on-set of vorticity in the condensate due to the cur¬ 
rent non-linearity. As in the standard case of a constant 
magnetic field, the effect of the density dependent gauge 
potential is to induce a rotation to the condensate, with 
the main difference that the onset of vorticity now de¬ 
pends on the number of particles in the condensate: the 
higher the density is, the more likely it is to have vortices. 

After introducing the mechanism by which an artificial 
gauge field acting on the condensate is generated, and 
deriving the most general expressions for the synthetic 
scalar and vector potentials, we introduce the conditions 
under which a constant magnetic field is generated. We 
derive the corresponding form of the generalised Gross- 
Pitaevskii equation, showing the appearance of a current 
nonlinearity. Finally we study analytically the conditions 
under which a vortex state appears as the ground state 
for the system in the Thomas-Fermi limit, comparing the 
result with the standard case of a rotating condensate. 


II. ATOMS IN ARTIFICIAL GAUGE FIELDS 

We consider a Bose gas of two-level atoms and 
we model the short-range interactions by a zero-range 
pseudo-potential. In addition, we include a coupling 
between the two levels, mediated by an external laser 
field. In the rotating wave approximation (RWA), which 
is valid when the frequency of the applied electric field is 
close to the atomic resonance, so that the terms in the 
Hamiltonian describing the atom-light interaction which 
are rapidly oscillating can be averaged to zero, in con¬ 
trast to the slowly oscillating ones [^, the microscopic 
Wbody Hamiltonian of the system is given by 

^ / 2 

n—1 ^ 

N 

+ ^'H\{n,e}, ( 1 ) 

n<l 

where the first term is the non-interacting Hamiltonian, 
with V (r) the external trapping potential, Pn the quan¬ 
tum momentum operator of the n-th particle, and m its 
mass. The identity matrices act on the sub¬ 
space excluding the particles n, Finally, the cou¬ 

pling t/a/(r) between the two atomic levels |1) and |2) is 
given by 



at the atom’s position r. These parameters, useful 
in order to make evident the unitary property of the op¬ 
erator, are related to the laser detuning A and the Rabi 
frequency k as A = F cos 0 and k = F sin 0. The detuning 
A is defined as the difference between the incident laser 
frequency and the frequency of the transition between 
the two atomic levels. The Rabi frequency k, character¬ 
izing the strength of the atom-laser coupling, is instead 
defined in terms of the matrix element di 2 = (1| d |2) of 
the atomic electric dipole operator between the two in¬ 
ternal states, and the amp litude E of the of the electric 
field, as K = di 2 • E/Ii [2^. 

The second term in the Hamiltonian m is the pairwise 
interaction between the particles, and has the diagonal 
form = diag [gn, gi 2 , gi 2 , 322 ] 5 (r„ - r^). The cou¬ 
pling constants are given by = Airh^aij /to, where Gij 
are the scattering lengths for the three different collision 
channels. 

In what follows we shall work in the weakly interact¬ 
ing limit, Pia^j 1 (with i,j = 1,2), and we consider 
the meanfield variational ansatz for the many-body wave 
function of the system: IF (ri, r2, ...r^r) = 
that is the symmetrized product of the single particle 
spinor wave function </>(r), satisfying the normalization 
condition J cPrcf>^(p = 1. We introduce the Lagrangian 
of the system, which reads as 


L = 


N 




{ihdt 




( 3 ) 


Substituting in eq.® the expression given above for the 
many-body wave function, we obtain the corresponding 
Lagrangian in terms of the single particle wave function 


Lmf= / a r 


{ihdt - Hmf) '0 


( 4 ) 


where 0(r) = '/Ncj){r) is the order parameter of the sys¬ 
tem (or wave function of the condensate), and we defined 
the single particle mean field Hamiltonian Hmf as 

HMF = ^®l + Uaf{v)+Vaa + V{v) (5) 

with I the 2x2 identity operator in the space of the 
atomic internal degrees of freedom. In eq.® Vaa de¬ 
scribes the mean field collisional effects, and is given by 

" 2 ( 0 

where 


_ fir / cosfi(r) e sinfi(r)\ , , 
2 sin fi(r) — cosfi(r) y 


= 5iiPi + gi2P2 (7) 

1^2 = gi2Pl + g22P2 ( 8 ) 


The light-matter interaction, eq. ®, is fully character¬ 
ized by three parameters, namely the generalized Rabi 
frequency F, the mixing angle 9{r) and the laser’s phase 


and Pi = \ipif' is the density of atoms in level |j), i = 
1 , 2 (with ipi the corresponding component of the order 
parameter). 
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In the weakly interacting limit, the coupling strength 
fir between the ground and the excited states is typi¬ 
cally much larger than the mean field energy shifts. This 
allows us to diagonalize the Hamiltonian in eq. m by 
treating the mean field interaction term as a small pertur¬ 
bation to the atom-field interaction. To order 0{iyd/fir), 
the eigenstates of eq. ([5]) take the form 

lx±) = lxi°^)±^lxi°^ (9) 

where Vd is the pertubative parameter, equal to 
10 0 

=-sin-cos-- 1 / 2 ) (10) 

and |x±^) are the unperturbed dressed states, having the 
form 

=cos^|l)-f e*bsin^| 2 ) 

IX-^) =sin^|l)-e*‘^cos^| 2 ) ( 11 ) 


The adiabatic approximation, used in order to derive the 
Hamiltonian in eq. m, is satisfied when the speeds of the 
particles are low enough that the time scale of their mo¬ 
tion is much longer than the time scale given by the Rabi 
frequency. Since the energy splitting between the inter¬ 
acting dressed state is proportional to the Rabi frequency 
in fact (and equal to fiT -|- (;/+ — i/_)), this condition pre¬ 
vents (or at least attenuates) the transition of the atoms 
from one state to the other. Upon substitution of eqs. 
(HH), (nni and m into ea. dTSl) . we obtain the following 
explicit forms for the vector and scalar potentials 

f% f 

A± = --(1 Tcos6»)V(;i-f ^p±V()) (16) 

2mW = ^(V6»)2-pfi2!^(V(/))2- 

^P+ cos6»(V(/))^ - 
HVe ■ V (^) (17) 

where we defined the coefficients f± as 


The condensate wave function \tj}{r,t)) can then be ex¬ 
panded in the interacting dressed state basis, eq. ®, 
as = Z)i={+,-} ^i(r,0lx*)- In order to cap¬ 

ture the dynamics of the ± component of the conden¬ 
sate we use the adiabatic assumption, according to which 
' 0 r|z(r, t) = 0 , and we consider the projection of the mean 
field Lagrangian in eq. (|3]) onto the subspace spanned by 
the corresponding (|x±)) dressed state. We obtain then 
the mean field Lagrangian for the relevant condensate 
component that, to order 0{i'd/fir) reads as: 


L± 



V'i {ifidt - H±)tjj± 


( 12 ) 


where 

+W±^ + i^± + V{r). (13) 

2m 2 

is the Hamiltonian describing the dynamics of the ± com¬ 
ponent of the condensate. Above we defined 

’^+ = 2 (2 +^2sm - 1 

If . 2 O 

zz_ = -(^.ism - + .2COS -j 

(14) 


while A± and W are effective vector and scalar po¬ 
tentials, respectively. Exploiting the identity relation 
I = |x+)(X-hI + it is easy to show that they 

get the form (see [ 1 ^fl^): 


A± = 
W = 


-(x±|p|x±) 

I(x-hIpIx-)I^ 


(15) 


/+ 


2 sin^ 6» 


\ , 

( ■ 2 O 

2 

COS -511 - 1 - gi 2 

2 

cos 


. 2^1 

-522 sin - 

/_ = -/+(1 O 2 ) 


(18) 


and with p± = I'i/'iP the density of atoms in the |x±) 
state. As we observe in eqs. (HU) and (II3, the vector 
and the scalar potentials also depend on the particle den¬ 
sity. In the following we analyse the important case in 
which, by properly engineered laser fields, the light-atom 
interaction induces an effective constant magnetic field, 
together with an extra nonlinear term that depends on 
the angular momentum of the condensate itself. 


III. PHYSICAL ASSUMPTIONS AND THE 
EQUATION OF MOTION 

We choose a particular laser configuration, such that 
the density dependent terms in the gauge potentials lead 
to an exotic nonlinearity in the Gross-Pitaevskii equation 
which is proportional to the angular momentum of the 
condensate and we will investigate the rotational prop¬ 
erties of such a system. With such setup, the lowest- 
order term in the vector potential corresponds to a con¬ 
stant magnetic field in the z-direction, whose effect on 
the physical properties of a condensate have been widely 
studied in the literature and are nowadays well under¬ 
stood 0 . It leads for instance to the nucleation of quan¬ 
tized vortices for sufficiently strong magnetic fields 
a situation that is equivalent to considering a rotating 
condensate in the rotating frame when the magnetic field 
qB = 2mflz, with Q the rotation frequency (and q = 1 
in our case). It is known in fact that once the angular 


2 m 
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velocity of the condensate (or the applied magnetic field) 
reaches a threshold value, a state with a vortex becomes 
energetically favourable with respect to a non-rotating 
one and, for strong magnetic fields, the Abrikosov trian¬ 
gular lattice of vortices is obtained. This type of ground 
state has been observed @ also in harmonically trapped 
clouds in the so called fast rotating regime, in which the 
centrifugal force nearly cancels the trapping potential. 
Based on these considerations, we will in the following 
address the question of how the density dependent gauge 
potential which, as said above, turns out to be propor¬ 
tional to the angular momentum of the condensate in our 
particular laser configuration, will affect the rotational 
properties of a condensate. 

There are several ways of creating such non-trivial ef¬ 
fective magnetic fields, and they typically depend on the 
details of the laser parameters. Here we choose to work 
in the symmetric gauge, where we obtain a constant mag¬ 
netic field by making the zeroth-order term of the vec¬ 
tor potential proportional to the radial distance from the 
origin of the reference frame, which we assume to co¬ 
incide with the minimum of any potential trapping the 
cloud. This can be achieved by properly designing the 
spatial dependence of the three parameters that are in¬ 
volved, namely A, n and the laser phase </>, in order to 
have —h/2 {1 ± cos9)V(j) oc r. To this aim we choose 
to work with the (-l-)-component of the condensate, and 
we choose a laser beam with a phase proportional to the 
polar angle ip as (j) = lip, where I is the quantum num¬ 
ber identifying the orbital angular momentum M carried 
by the electromagnetic field. It can be shown that, by 
making the assumption k <C A, the sought spatial de¬ 
pendence for the zeroth-order term is obtained, once we 
tune K, A in such a way that k/A ^ r. It is worth noting 
that if we had chosen to work with the (—)-component 
of the condensate, the dynamics of the system would be 
the same provided we change the sign of the detuning, 
and we exchange gn O (722 • 

The condition described above can be realised by de¬ 
signing the laser beam in such a way that the Rabi fre¬ 
quency is linear in r in the plane in which we confine the 
cloud, and equal to K(r) = kqt where kq is a constant and 
also assuming that the detuning is constant in the same 
plane. These considerations are clear once we substitute 
the following relations in eq. (HZl: 


sin0 ~ k/A 
cos9 ~ 1 — k^/2A^ 

9 ~ k/A (19) 


that hold true in the limit k A. In this limit, the 
perturbative condition <C HT reduces to Vd ^ hA, 


and the expressions for the potentials become 


A - 


712 
A2 A 


pV(j) 


( 20 ) 


2,„W -(v£) 

-'‘X (’i) ■[("'!) "+1^ 

together with 


~ -A 1 + 


2A2 J 


511 « 

- ^P-712^P 


( 22 ) 

(23) 


in which we defined 712 = {gn — 512 ) /4, and we dropped 
for brevity the (-I-) subscript, since from now on we are 
dealing just with the atoms in their |x+) internal eigen¬ 
state. 

To justify the assumption we made for the ratio k/A oc 
r, and in particular about how it is restricted by the 
condition k <C A, we dehne a ^ 1 such that k/A = 
kot/A < l/a. This condition is fulfilled in the region 
r < A/Koa = R that also defines the region of validity of 
our model, and can be modified by properly tuning the 
constants A and kq. 

In order to arrive at a Gross-Pitaevskii equation for 
the order parameter of the condensate, we minimize the 
action 


S= dtL+. 


(24) 


with respect to ?/>*. This leads to the equation of motion 


..d,p 

'’"m = 


(P “ A)^ hn 

+ V {r) + W + Y +'^+ 


2m 


-b 


(ar.jr+a2.j2) + —(v -) 

+ ^P + 7i2p(£) V' (25) 

in which the nonlinear current ji and the j 2 term defined 
below appear, along with the coupling vector fields ai 
and a 2 , with the form 


ai 

ji 

a2 

h 


712 _ 

A2 A 

h 

2im 
n 712 

A A 

— VV'+ lAVi/-*] . 
2im 


iIj — i- 

(-z) 


tp* -t/j* [V+ i 


(26) 
7) 
(28) 
(29) 


Using the explicit expressions, and (p = 19, the 
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Gross-Pitaevskii equation reads as 


'’‘m = 


M / Ko \ ^ r T r / \ 

+ 5 ;; (a) 


2m 




(I) 


- ( gii + 2712^ I 




(30) 


where L 2 = [V'^ ~ ^Sy) orbital angular mo¬ 

mentum operator. As expected, the linear term in ea. dSOl) 
is identical to the Hamiltonian describing the motion of 
a charged particle in a constant magnetic field (or of a 
rotating condensate). The relation between the effective 
magnetic field B, the angular velocity fl and the param¬ 
eter of our system, is now given by 

B = 2 mn = y . (31) 

The first bracket in eq. ( 1501 ) describes the dynamics of 
a particle, subjected to the effective magnetic field eiD 
whose intensity is proportional to the orbital angular mo¬ 
mentum M carried by the laser beam, and confined by the 
potential V (r), while the second bracket is the resulting 
nonlinear term of first order in the atom-atom interac¬ 
tion. In the latter, the second term contains the stan¬ 
dard mean field term proportional to the density 511/9. 
We will drop the last term in what follows, since it is of 
order 0(k/A)^, and corresponds to a scalar interaction 
much smaller than 5 in the centre of the trap. The goal 
is to study the role of the nonlinear term proportional 
to the orbital angular momentum Lz, since the remain¬ 
ing terms only modify the effective mean field coupling 
constant. In order to isolate the effect of the density de¬ 
pendent gauge potential we envisage a situation where 
we counteract the linear zeroth order gauge potential by 
rotating the system at the frequency H = — ^ (^) , 
along with the laser apparatus. This might be exper¬ 
imentally challenging |30|, but it allows us to focus on 
the effect of the density dercndent gauge potential, from 
which novel phenomena |25| is expected to arise. The 
resulting Gross-Pitaevskii equation, whose study will be 
the object of the next section, takes finally the form 




CpLz -f V (r) -|- gp 




(32) 


where we have defined 


IV. THE VORTEX STATE 

In the previous section we derived the Gross-Pitaevskii 
equation for a condensate with an angular momentum 
non-linearity. We will next focus on the physical con¬ 
ditions for the nucleation of vortices in the condensate. 


In order to extract any new phenomena deriving from 
the nonlinearity, we will consider the case of a state con¬ 
taining a single vortex state carrying one quantum of 
circulation. 

An interesting point that is worth emphasizing is that 
there are two different ways in which we can look at the 
nonlinear angular momentum term. On one hand, it can 
be seen as a modification of the mean field coupling con¬ 
stant, which depends on the local angular momentum of 
the condensate. This results in a breaking of the rota¬ 
tional invariance of the system, enhancing or weakening 
the effective interaction felt by the atoms, depending on 
the direction of rotation of the condensate. On the other 
hand, the nonlinearity can also be understood as a den¬ 
sity modulated angular velocity of the cloud. In order to 
see to what extent the rotational properties of the con¬ 
densate are modified, we start the discussion considering 
this latter interpretation, looking for the critical value of 
the C parameter in eq. dSSD, that makes a vortex state 
of the condensate energetically favourable with respect 
to its non-rotating solution. To this aim we consider a 
BEG in a highly anisotropic trap, such that its dynam¬ 
ics in the transversal direction is frozen and the system 
is quasi-two-dimensional, and confined in its plane by an 
isotropic harmonic potential of frequency ujq. We assume 
the condensate uniform in the z direction, with thickness 
Z. The original order parameter can then be rescaled 
as il}{v)l'/Z^ where '(/'(i’) is now two-dimensional, and 
normalized in such a way that f = A^- All the 

previous expressions for the Hamiltonian JI± and the po¬ 
tentials A±, W hold true, with the meanfield coupling 
constants and the C parameter rescaled as 5 —>■ g/Z and 
C —>■ C/Z. We consider the condensate composed by 
a number of atoms sufficiently large so that its physi¬ 
cal properties are well described working in the Thomas- 
Fermi limit [3l|. 

In this regime, the width R of the condensate is much 
larger than the healing length ^ so that, because of the 
large r behaviour of the azimuthal velocity field, the extra 
energy associated with the presence of a vortex (com¬ 
pared to the energy of the same non-rotating system) 
can be estimated by considering just the corresponding 
kinetic term [m giving Ev = 7rp(0)^log(0.888i?/^), 
where p{Q) is the density of the non-rotating condensate 
at the center of the trap. In the standard situation in 
which only the linear rotational term —ilLzip is in the 
Gross-Pitaevskii equation, the vortex state becomes en¬ 
ergetically favourable when the energy Ey = Ey — ^{Lz) 
of the excitation in the reference frame co-rotating with 
the condensate is less than the energy corresponding to 
the non rotating cloud, with (Lz) the value of the angu¬ 
lar momentum in the single vortex state. This condition 
leads to the critical angular velocity 

7 ^ ofc 

^ (Z^ ^ (0-888i?/^) (33) 

where (Lz) = / cPrip*Lz'ip = ^p{0)TrhR^, having used 
the expression p{r) = p{0) (l — r'^/R'^) for the Thomas- 
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Fermi density profile. Motivated by this, we note 
that, with the nonlinear term —Cp{v)Lz'4! in the Gross- 
Pitaevskii equation, we can estimate the critical value 
Ccr in the same fashion, obtaining: 


Ccr 


3 ricr 
2 ^ 


(34) 


having used {pL^) = f (Prp{v)'ip*L= ^hp^{Q)E?. 
Apart from a constant factor then, the critical value of C 
is given by Ccr ^ same result can be obtained 

by looking at the new nonlinearity as a modification of 
the mean field coupling constant. This effect is clear once 
we substitute the ansatz ^(r) = /(r)e*‘^ for the order pa¬ 
rameter into eg. (1321) . with the characteristics of a vortex 
state with angular momentum h per particle. In this way, 
we get the GP equation for the modulus f{r) of the wave 
function, with the mean field coupling constant changing 
a.s g ^ g' = g — hC. The effective coupling constant 
then is stronger or weaker depending on the combina¬ 
tions of the scattering length differences and the orbital 
angular momentum of the incident laser beam, which can 
contribute both with a positive or negative sign. 

We then see that the nucleation of a vortex changes 
abruptly the properties of the condensate since it changes 
the mean field interaction strength, which also affects 
physical quantities such as the healing length and there¬ 
fore also the shape and size of the vortex core. Moreover, 
the modification of such interaction strength shifts the 
energy of the system at hand, for which the energy of 
the non-rotating state can be taken as a reference value, 
and in respect to which all the excited states of the con¬ 
densate are defined. 

In the case g' < g, the abrupt drop of this back¬ 
ground energy could make the vortex solution energet¬ 
ically favourable compared to the non rotating one. This 
happens when this energy shift is larger than the energy 
Ev of the vortex, see Fig. 1(a). 

In the TF limit, the total number of particles is related 
to the chemical potential by the equation 


N = 




Uo mwg 


(35) 


having indicated with Uq the generic value of the mean 
field coupling constant. Solving eq. (l35ll for p, we get 


p = 


itiujqUo 

TT 


1/2 

NC2 


(36) 


and since by definition p = ^ we get, by using eg. (1551) 

cN 2 

Eo{N) = / dN'piN^) = -pN. (37) 

Jo o 

The energy difference AEq between the two non-rotating 
states with values Uq = Uq — AUq and Uq for the mean 
field coupling constant, is equal to: 

AEq = Eq(Uq) - Eq(Uq - AUq) = ^^°}^°^ AUq (38) 

oUq 


(a) (b) (c) 

C> Qr C = Qr C < Qr 


E 



FIG. 1: Illustration of the conditions for the onset of vorticity 
in the condensate, when the energy of the vortex state is low¬ 
ered by the effect of the laser field incident on the condensate 
itself. The solid long line represents the reference non rotating 
state energy and the solid short one the energy of the vortex 
state, (a) When the redefinition of the meanfield coupling 
constant g' ^ g — tiC, due to the presence of the vortex in 
the superfiuid, induces a shift AEq of the non-rotating state 
energy Eq , which is larger than the energy Ev of the vortex it¬ 
self, then the rotating state becomes energetically favourable 
compared to the non rotating one. (b) The critical condition 
happens when Ev equates the energy shift of the non rotat¬ 
ing state. In the case (c), the energy of the vortex state is 
lowered, but not enough to induce its nucleation. 


for AUq <C Uq. Differentiating eg. (1551) with respect to 
Uq, we have dp/dUo = l/2p/Uo so that: 

AEq = (39) 

This energy shift has to be compared with the en¬ 
ergy of the vortex state Ey (Uq — AUq) = Ey (Uq) — 
[dEy(UQ)/dUQ],AUQ. Since p{0) = p/Uq, we get 
dpiO)/dUQ = -\l2plUl, with \dEy(UQ)ldUQ\ AUq = 
— 1/2 (p/Uq) (AUq/Uq) being a higher order correction to 
the vortex energy Ey (Uq), and can be neglected. Equat¬ 
ing Ey (g) with the value AEq calculated for AUq = HC, 
that in ordinary cold atoms experiments fulfills the con¬ 
dition hC <C g, the critical value for the C parameter 
takes the form 

C„ = 3^ (40) 

that is the same result as in eq. ([551). up to a multiplica¬ 
tive factor of 2. This discrepancy can be ascribed to the 
fact that the result in eg. (|55)) represents an estimation 
for Ccr, since it is obtained by making an equivalence 
between the system at hand and a rotating condensate, 
that actually are two different physical situations. 

This result tells us that the effect of the synthetic gauge 
potential is in fact to induce a rotation in the condensate, 
self-boosted by the total number of atoms composing the 
condensate itself. The higher this number is, the easier it 
is to induce the nucleation of vortices and so to put the 
cloud into rotation. Remembering the definition of the 
parameter C in terms of the physical quantities charac¬ 
terizing the system, we obtain the critical value 

/^-(O) = (-]' log (0.888R/O (41) 

«712-R \kq J 
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of the particles number for which a vortex state is 
favoured. From the expression in eq. dnj we see that 
there is significant flexibility in reaching the required crit¬ 
ical combination of parameter values, due to the the fact 
that the orbital angular momentum of the light beam t, 
the scattering length difference 712, and also to some ex¬ 
tent the density of the cloud, can all be relatively easily 
changed in an experiment. For example, we can choose 
the scattering length difference an — 012 to be 100 nm 
by using Feshbach resonances if necessary, a Rabi fre¬ 
quency of the order of 10“^ Hz at the border of the cloud, 
and a value for the homogeneous laser detuning from 
the atomic resonance such that kqR/A = 0.1, satisfying 
the assumption k <C A throughout the condensate (with 
a = 10 at the border of the cloud). A metastable Helium 
gas could be particularly useful in our case because of the 
small mass of the particles 10“^^kg). Considering an 
angular momentum carried by the laser beam character¬ 
ized by the quantum number I ~ 100 ( 3 ^ . we obtain a 
critical value for the particle density at the centre of the 
trap Pcr(O) of the order of 10^^ atoms/cm^, for which the 
vortex state is energetically favourable compared to the 
non rotating one. This density can be reached in current 
experiments for which the condition of weakly interacting 
atoms pc? 1 is still fulfilled. 

The result in eq. m, obtained for the case of a vor¬ 
tex state carrying a single quantum of circulation, can 
be easily generalised to the case of a vortex with higher 
angular momentum. Similar results are also expected for 
states with more than one vortex. 


V. CONCLUSIONS 

In this paper we have studied the energetic properties 
of a singly quantized vortex state in a Bose-Einstein con¬ 


densate which is subject to a density dependent gauge po¬ 
tential. Nontrivial nonlinear dynamics emerges through 
the collisionally induced detuning for the light-matter 
coupling, which manifests itself as a nonlinearity where 
the density times the angular momentum of the conden¬ 
sate appears. We have chosen to study one particular 
setup using an incident laser beam with an orbital angu¬ 
lar momentum, which couples the two internal levels of 
the atoms. This is only one possibility; there are many 
other setups which can be used which do not explicitly 
rely on a symmetric gauge ©[II but will provide a 
uniform magnetic field for the zeroth order term. Nev¬ 
ertheless, the main conclusion still remains, the onset of 
a rotating state via the effective magnetic field, will de¬ 
pend on the density of the condensate in combination 
with other parameters, in stark contrast to the standard 
situation of vortex creation. Apart from the fundamen¬ 
tal interest in rotating superfluids, the results presented 
here also shine light on the unprecedented flexibility in 
creating novel types of nonlinear dynamics for the Bose- 
Einstein condensate. The current non-linearity used in 
this work opens up many new questions regarding the 
superfluid dynamics. Eor instance, drag forces and the 
vortex dynamics including vortex interactions may well 
be dramatically altered due to the current non-linearity. 
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